tt = 0:dt:Dt;
mtt = 0:dt:mDt;
ptt = 0:dt:pDt;

figure
plot(tt,(R(1:t+1)-RL)/1000,'r');
hold on 
plot(mtt,(mR(1:mt+1)-RL)/1000,'g');
plot(ptt,(pR(1:pt+1)-RL)/1000,'b');
xlabel('t/s');
ylabel('h/km');

figure
plot(tt,U(1:t+1),'r');
hold on 
plot(mtt,mU(1:mt+1),'g');
plot(ptt,pU(1:pt+1),'b');
xlabel('t/s');
ylabel('u m/s');

figure
plot(tt,V(1:t+1),'r');
hold on 
plot(mtt,mV(1:mt+1),'g');
plot(ptt,pV(1:pt+1),'b');
xlabel('t/s');
ylabel('v m/s');

figure
plot(tt,W(1:t+1),'r');
hold on 
plot(mtt,mW(1:mt+1),'g');
plot(ptt,pW(1:pt+1),'b');
xlabel('t/s');
ylabel('w m/s');

figure
plot(tt,ALPHA(1:t+1),'r');
hold on 
plot(mtt,mALPHA(1:mt+1),'g');
plot(ptt,pALPHA(1:pt+1),'b');
xlabel('t/s');
ylabel('alpha rad');

figure
plot(tt,BETA(1:t+1),'r');
hold on 
plot(mtt,mBETA(1:mt+1),'g');
plot(ptt,pBETA(1:pt+1),'b');
xlabel('t/s');
ylabel('beta rad');